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Abstract A brief review is first presented of attempts to predict stable niultiquark states within 
current models of hadron spectroscopy. Then a model combining flip-flop and connected Steiner trees 
is introduced and shown to lead to stable multiquarks, in particular for some configurations involving 
several heavy quarks and bearing exotic quantum numbers. 
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1 Introduction 



The problem of the existence of exotic hadrons is even older than the quark model. In the early 60s, 
the speculations on a possible Z baryon concerned states with unusual strangeness, as compared to 
the hyperons. Nowadays, genuine exotics are defined as hadrons whose quantum numbers cannot be 
matched by ordinary quark-antiquark (qq) nor three-quark (ciqq) configurations of the simple con- 
stituent picture. Gluonium, or hybrid, or multiquark states having the same quantum numbers as (gg) 
or {qqq) are sometimes referred to as "cryptoexotics". 

There are at least two major and recurrent problems in this field, which might often discourage 
the newcomers: the ups and downs of the experimental searches, and, for theorists, naive enthusiasm 
alternating with excessive scepticism. 

The prospects are, however, encouraging. The recent experiments at Belle, BaBar, Cleo, Fermilab, 
etc., have shown that new heavy hadrons can be identified, such as the ric{2S), the f?^ or the Ar(3872). 
On the theory side, the approaches based on lattice QCD, QCD sum rules, and even AdS/QCD usefully 
supplement the studies based on constituent models. 

This survey is organised as follows. In Sec. [2] the experimental situation is summarised. In Sec. [3] 
the early and recent speculations about multiquarks are reviewed. After a discussion about the chromo- 
electric model in Sec. |4] the Steiner-tree model of confinement is described in Sec. [5] and the results 
based on this picture are displayed in [6] before some conclusions in Sec.[7j 
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2 Some experimental results 

The Z baryon was tentatively predicted by analysing the data on KN scattering, and never firmly 
seen. For a review, see early issues of "Review of Particle Physics", such as [1] or others as listed in the 
latest one |38| . 

Then came the baryonium. See, e.g., f36]. There were puzzling peaks in the p cross-sections, and 
also in the inclusive 7 spectrum of — > 7 + X at rest. The most intriguing indication came from an 
experiment by French et al. fT9] in which a narrow state of mass about 3 GeV was seen decaying often 
through another baryonium of mass 2.2 GeV. These claims provided the main motivation to build the 
low-energy LEAR facility, as a side-project of the ambitious SppS collider project at CERN. None of 
the baryonium candidates were ever confirmed at LEAR and elsewhere. However, some enhancements 
have been observed recently in the mass spectrum as studied in B decay at Belle and BaBar [15], 
or in J/V' decay at BES |TS] . 

More recent, and perhaps not yet fully settled, is the issue of possible light pentaquarks, as suggested 
in some speculative developments of chiral dynamics. The experiment by Nakano [35] is commented 
upon in some other contributions to this meeting. It had at least the merit to look at new spectroscopy 
against the fashion of that time. More surprising, as discussed, e.g., in |55i i51, is the wave of followers: 
data on tape were hastily analysed, that nobody had ever the curiosity to look at; lattice and QCD 
sum rules calculations were quickly published, with uncertainties even on the basic quantum numbers 
and no clear distinction between genuine bound states or resonances and states artificially quantised 
by finite-volume effects. Fortunately, there are nowadays more reliable lattice or sum-rules studies, as 
reported in this conference, or, e.g., in [111 HO] . 

The case of dibaryons is more dilute in time. There are persisting claims, which have neither been 
firmly confirmed nor definitely ruled out. Examples are |52t I31| . There were detailed predictions in 
the late 70s (see, e.g., ||37j), using models elaborated for baryonium, but most states were due to an 
artificial bag confinement neglecting the mere "fall-apart" decay into two baryons |28| . 

3 Early and recent mechanisms 

3.1 Duality 

We refer to the review by Roy [17]. Duality was devised to give consistency to s-channel and i-channel 
pictures of hadronic reactions. In ttTV scattering, for instance, summing over all baryon resonances 
should be equivalent to the cumulated meson exchanges in the crossed channel. Rosner |45l pointed 
out that the crossed-channel partner of meson exchange in NN scattering consists of mesons made of 
two quarks and two antiquarks and preferentially coupled to NN: the baryonium was invented! Note 
already the Pandora-box syndrome: once baryonium is accepted, duality applied to baryonium-baryon 
scattering implies the pentaquark |i7] . 

3.2 Colour chemistry 

The name is due by Chan H.M. and his collaborators [13]. It describes the art of using clusters to 
simplify the quark dynamics and extrapolate toward higher configurations. There are, however, several 
levels: 

1. The diquark was invented to simplify the picture of baryons. As compared to the symmetric 
quark model, the quark-diquark spectrum does not include the configurations in which both degrees 
of freedom are excited, such as a; x y exp(— a(a;^ -I- y^)/2) in the specific harmonic oscillator model, 
with X and y being the Jacobi variables. See, e.g., [?]• 

2. This diquark with colour 3 has been used to build new mesons, as diquark-antidiquark states [27] . 
A warning was issued, however, that taking too seriously the diquark could lead to unwanted " demon" 
deuterons |21lM This holds for recent extensions to the heavy-quark sector [331, in which no caution 
is cast that a \cq) — (cq) picture of some new mesons could lead to unwanted [{cu){cd){cs)] dibaryons 
below the (ccc) -I- (uds) threshold. 



^ There is an error on the quantum numbers of the {qq)^ state, but the issue remains. 
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3. A further speculation consists of assigning these clusters to some higher representation of the 
colour group. In [M], colour-sextet diquarks were invoked to predict massive mesons whose both 
mesonic and baryonic (NN) decays are suppressed. This was also proposed for the late pentaquark |49| . 

3.3 Chromomagnetism 

In the mid-70s, the chromomagnetic interaction jl7|. 

Vss^-AY.'^Xr^.Xrcr..^,, (1) 

i<j 

inspired by the Breit-Fermi interaction in atoms, was shown to account for the observed hyperfine 
splittings of ordinary hadrons, including the S — A mass difference. There is a relativistic analogue 
in the bag model. In |26| . JafFe pointed out the H{uuddss) in its ground state has a chromomagnetic 
interaction whose cumulated (attractive) strength is larger than for the lowest threshold A{uds) + 
A{uds), due to the non-trivial colour-spin algebra in ([T|) 

(Ar'.Af)...a,>^ = 2(...>„,. (2) 

If the quark masses m,; are identical (SU(3)f symmetry) and if the short-range correlation factors 
{5^'^\vij)) are assumed to be the same in multiquarks as in ordinary baryons. then the H is predicted 
about 150 MeV below the threshold. The H was searched for in more than 15 experiments, with mostly 
negative results, including double- yl hypernuclei. An enhancement above the AA threshold cannot be 
excluded jH]. 

In 1987, Lipkin, and independently Gignoux et al. [351 '^^^'SS^ pointed out that the same mechanism 
would bind an exotic pentaquark (the word was first used in this context) 

(Qqqqq) < (Qq) + (qqq) , (3) 

where Q = c, b and (qqqq) is a combination of u, d, s, in a SU(3)t flavour triplet. The same binding of 
— 150 MeV is found if Q becomes infinitely heavy. A search for the 1987- vintage of the pentaquark in 
an experiment at Fermilab turned out inconclusive [3J|3]. Other candidates have been listed [50 . 

Further studies of the chromomagnetic model |46l [501 EUl Rt8j indicated that the corrections tend 
to moderate and even to cancel the binding effect due to the group-theoretical result (|2|, namely 

— There is more kinetic energy in a dibaryon than in two baryons (5 Jacobi variables vs. 4), 

— SU(3)f breaking is not favourable. The A resists much better than H. 

— The matrix elements {S^^'^{rij)) are significantly smaller in multiquarks than in ordinary hadrons. 

Nevertheless, the simple chromomagnetic Hamiltonian ([ij might still reveal some surprises. For 
instance, if diagonalised with realistic values for the flavour dependent strength factors {S{rij) /{mirrij), 
it provides the (ccqq) lowest configuration with J^'-^ = l"*""*" the main features required to describe the 
A:(3872) US [12]. 

4 Symmetry breaking in the naive chromoelectric model 

If the chromomagnetic mechanism turns out disappointing, it is natural to consider the chromoelectric 
interaction. Its major property, besides confinement, is flavour independence, as in QCD, giuons couple 
to the colour of quarks, not to their fiavour. The situation there is reminiscent from the physics of 
exotic atoms, where e~, /i~ and p create the same electric field. 

There is an abundant literature on the stability of few-charge systems such as {m^ , , , ) 
as a function of the constituent masses rUi. See, e.g., |5] and refs. therein. In particular, it has been 
proved that the positronium molecule Ps2 is stable in the limit where annihilation is neglected. This cor- 
responds to equal masses. Then breaking charge conjugation leads to a more stable {M~^ , M+, m~ , m~) 
whose prototype is the hydrogen molecule. However, breaking particle identity gives (M+ , m+ , AI" , mT ) 
which becomes unstable against dissociation into (M+,M~) -I- [rri^ ,m~) if mjM < 2.2 (at most a 
metastability with respect to {M^,m~) + {m^,M~) can be envisaged). 
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But the Coulomb character of the interaction does not matter, and the trend can be extended 
to quark models with flavour independence. For instance, it was noticed rather early that {QQqq) 
becomes bound below the (Qq) + {Qq) threshold if the mass ratio becomes large enough [51 "M^ . For a 
discussion of more recent studies, see |53j . 

It remains to understand why, contrary to the positronium molecule with Coulomb forces, the equal- 
mass {QQQQ) system is not found to be stable in simple quark models based on the colour-additive 
ansatz 

^--^E^^^^-?^ (4) 

which is normalised so that v{r) is the quarkonium potential. Let us consider a class of Hamiltonians 
supporting at least one bound state, and submitted to the constraint that they have the same average 
strength g 

H[g]=Yl P*'2m + 9^An,) , = ^^^^^^5 • (5) 

i i<j 

Using the variational principle, it is easily seen that the symmetric case [g]s where gij = g gives 
the highest energy, and that the most asymmetric sets of coupling lead to the lowest energie^ Let us 
now compare in table ^ the trend of couplings for two isolated mesons or atoms, in the positronium 
molecule and in a tetraquark governed by Q. The spread is measure from the variance, and both 3 — 3 
and 6 — 6 colour configurations are considered. 



Table 1 Spread of couplings for the positronium molecule and the tetraquarks in the simple colour-additive 
model, and in their thresholds. 



(abed) 


v{r) 


3»j 
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(l,3)+(2,4) 


— 1/r, r 


{0,0,1,0,1,0} 


1/3 


0.22 


PS2 


-1/r 


{-1,-1,1,1,1,1,1} 


1/3 


0.89 




— l/r,r 


{1/2,1/2,1/4,1/4,1/4,1/4} 


1/3 


0.01 


[{qq)6{qq)(i] 


— 1/r, r 


{-1/4,-1/4,5/8,5/8,5/8,5/8} 


1/3 


0.17 



The molecule Ps2 is stable because is uses more asymmetry than two isolated atoms. This is not 
exactly the way it is taught in textbooks on quantum chemistry, but this property is equivalent to the 
statement that in Ps2 the two atoms become enough polarised to ensure their binding. On the other 
hand, the balance is unfavourable for both [{qq)3{qq)3] and [{qq)eiqq)6\- 

This means that multiquark stability is penalised by the non-Abelian character of the colour algebra 
with restricts the variation of the strength of the potential when colour is modified, unlike the Abelian 
case where it can flip sign. 

In explicit quark model calculations, such as |29l I53|. short-range and spin-dependent terms are 
added, but there is no stable tetraquark in the case of equal masses, if this colour-additive rule is 
applied. It is found that is some models, (ccud) is marginally stable. 



5 Flip-flop and Steiner-tree confinement 

It has been anticipated by many colleagues that the confining potential needs not to be pairwise. Artru, 
Dosch |5J [TH] and many others have suggested that if the quark-antiquark potential of quarkonium 
reads v{r) = err at large distances, the analogue for baryons is the so-called "Y-shape" potential (see, 
Fig.0 

3 

Vy = (J min ^ ■ (6) 

s ^ — ^ 

i=l 



^ To be rigorous, one can compare only sets of coupling [g] that are aligned with the symmetric case [g]s in 
the space of coupling constants: then the farther of [g]s, the lower the energy. 
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This potential is found, e.g., in the adiabatic version of the bag model [23 , suited for heavy quarks, 
in which confinement is enforced only for gluons, and in various variants of the flux tube model. 
Estimating Vy is nothing but the celebrated problem of Fermat and Torricelli, which, in turn, is linked 
to a theorem by Napoleon. With the notation in Fig. [T] the length of the baryon string is given by 
VV/ct = Y2i IN'f^ill = Ika^all) is equal to the distance of a quark, say W3, to the auxiliary point -Ug 
completing an external equilateral triangle with the two other quarks, vi and V2, and often referred 
to as Melznak's point in the literature on minimal paths. Napoleon's theorem states that if the three 
external equilateral triangles are drawn, V2V3v[ and tisWiWjjtheir centres form an equilateral 

triangle, an interesting example of symmetry restoration if the initial triangle V1V2V3 is asymmetric. 



• O 



Fig. 1 Quark-antiquark (left) and three-quark confinement in the minimal string limit. 

Unfortunately, estimating the baryon spectrum using this Y-shape potential |44| gives results that 
are very similar to these from the colour-additive model, which reads for baryons (it is known as the 
"1/2" rule) 

2 (^12 + 7-23 + T'ai) • (7) 

To extend this potential to the case of niultiquarks, the recipe was first a little empirical, but it has 
been later endorsed by lattice QCD [41 and AdS/QCD [6 . For tetraquarks, it reads 

Vi = min [cr(ri3 + r24), cr{ri4 + r23), Vyv] , (8) 

a minimum of the "fiip-fiop" term (most economical double quark-antiquark pairing) and a connected 
Steiner tree linking the quark sector to the antiquark one and generalising the Fermat-Torricelli prob- 
lem beyond three points. In the case of planar tetraquarks, the length of the Steiner tree can be 




Fig. 2 Generalisation of the linear potential of mesons to tetraquarks: the minimum is taken of the flip-flop 
(left) and Steiner tree (right) configurations. 

estimated by iterating Napoleons's construction, as shown in Fig. [3) In the case where the quarks and 
the antiquarks do not belong to the same plane, the Melznak point W12 lies somewhere on a circle C12 
whose axis is the line V1V2 joining the quarks, and similarly for W34 on the circle C34 in the antiquark 
sector. To estimate the connected potential Vyy one can replace the minimisation over the location of 
the Steiner points si and S2 in ||sii'i|| -I- ||sii'2|| + ||'S2f3|| -I- ||s2W4|| by a maximisation of ||wi2U'34|| over 
the Melznak points, i.e., the maximal distance between the two circles. See Fig. [3] right. This leads 
to a bound on Vyy, that can be shown to hold for the overall interaction V4, and reads (mi2 is the 
middle of the quarks, and 77134 of the antiquarks) 

/3 

Vi/a < —i\\viV2\\ + \\v3Vi\\) + II 777 1 277734 II , (9) 

This very crude inequality enables one to bound the 4-body Hamiltonian in a form that splits into 
three independent pieces, one per Jacobi coordinate, and to demonstrate the stability of the tetraquark 
in some hmiting mass configurations [lOj . 




W3i 




Wl2 



Fig. 3 In the planar case (left), the length of the Steiner tree for planar tetraquarks, Wy/o" ~ ||siWi|| + 
||sif2|| + IIS2W3II + ||s2U4||, is equal to the distance between the two Melznak points «;i2 and U134. In space 
(right), the length of the minimal Steiner tree is is the maximal distance of the circles C12 and C34. Every point 
of C12 makes an equilateral triangle with the two quarks, and similarly for C34 for the antiquarks. 



This tetraquark potential, that combines flip-flop to a connected Steiner tree, can be easily extended 
to higher conflgurations. For the pentaquark, the flip-flop interaction can be read on Fig. [l] if a 
permutation over the quarks is implied. The connected diagram is shown on Fig. 4] (left), again with a 
permutation over the quarks [43j . It looks like an antibaryon, with two antiquarks replaced by a pair 
of quaxks in a colour 3 state. 





Fig. 4 Left: connected contribution to the pentaquark confining potential. Centre and right: contributions to 
the dibaryon potential 

For the dibaryon, the minimum should be taken of two separated Y trees and a connected tree, 
with all permutations, as shown in Fig. |4] (centre and right). 

For three quarks and three antiquarks, the choice for the minimum is between three mesons, a 
meson and a tetraquark, a baryon and an antibaryon, or a connected diagram. SeeFig. [5] 



> 





Fig. 5 Contributions to the potential of three quarks and three antiquarks. 



6 Results 



Needless to say that calculating the 4-, 5- or 6- body problem with this potential is more delicate 
than with ordinary pairwise models. This can be done with either the hyperspherical formalism or an 
expansion into a basis of correlated Gaussian functions, and in both cases, the matrix elements have 
to be computed numerically for each call of the potential energy. 

A first attempt concluded to the "absence of exotics" within this model [13 , but the 4-body wave 
function was too crude, and the study was restricted to the case of equal masses. As already mentioned, 
it can be demonstrated analytically that the tetraquark {QQqq) is bound within this potential in the 
limit of very large quark-to-antiquark mass ratio M/m [TUJ. A numerical estimate indicates binding for 
any value of the mass ratio M/m, including M = m, while for (QqQq), binding is quickly lost [54^ . The 
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pentaquark (qqqqq) is also bound, either for equal masses or in the limit where the antiquark or one of 
the quarks becomes infinitely heavy j43]. Javier Vijande, Alfredo Valcarce and I are presently looking 
at the 6-body problem, either the hexaquark (QQQqqq) or the baryon-antibaryon one, (QQQqqq), 
and we have indication of binding for several values of the mass ratio M/m. 

Note, however, that in the above studies, no account is taken of the role of short-range forces, 
such as the Coulomb-like interaction in the central potential or the spin-dependent forces. Also, when 
equal masses are involved, it is understood that the quarks remain distinguishable. In this Born- 
Oppenheimer like scenario, for each position of the quarks and antiquarks, the lowest configuration of 
the gluon fiux is selected, but the various string topologies correspond to different colour couplings and 
thus to different symmetry patterns for the spin and orbital wave functions. It is hoped nevertheless 
that more realistic estimates could be elaborated on the basis of this model of confinement, with full 
account for the Pauli principle, short-range forces, spin-dependent corrections, etc. 



7 Conclusions 

The lowest-order three gluon diagram is suppressed in the baryon dynamics, due to its colour structure, 
though it might play a role in the spin-dependent quark distributions at very high momentum, as 
stressed recently by Mitra [35', one of the founders of our few-body community. 

The non-Abelian nature of the dynamics is also crucial at large distances, with presumably the origin 
of confinement, and the multi-body character of the confining interaction. For baryon spectroscopy, the 
Y-shape interaction does not change very much the predictions, as compared to empirical models based 
on pairwise potentials. However, the new multi-body interaction, which is now supported by lattice 
QCD, give more attraction for multiquarks, and leads to stable configurations. This should stimulate 
new experimental investigations, in particular in the heavy quark sector. In particular, the search for 
(QQqq) tetraquarks could be combined with the search of (QQq) baryons with two heavy quarks, with 
the same triggers and analysis devises. 

Acknowledgements It is a pleasure to thank the organisers of this Conference for the opportunity of many 
stimulating discussions in a beautiful framework. 
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